PfafRan and fragmented states at i/ = | in quantum Hall droplets 
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When a gas of electrons is confined to two dimensions, application of a strong magnetic field 
may lead to startling phenomena such as emergence of electron pairing. According to a theory this 
manifests itself as appearance of the fractional quantum Hall effect with a quantized conductivity at 
an unusual half-integer «^ = f Landau level filling. Here we show that similar electron pairing may 
occur in quantum dots where the gas of electrons is trapped by external electric potentials into small 
quantum Hall droplets. However, we also find theoretical and experimental evidence that, depending 
on the shape of the external potential, the paired electron state can break down, which leads to a 
fragmentation of charge and spin densities into incompressible domains. The fragmentation of the 
quantum Hall states could be an issue in the proposed experiments that aim to probe for non-abelian 
quasi-particle characteristics of the v ~ j quantum Hall state. 



PACS numbers: 73.21. La, 73.43.-f, 71.10.Pm, 85.35.Be 
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I. INTRODUCTION 



The discovery of the quantum Hall (QH) effect at Lan- 
dau level filling factor ^' = | in the two-dimensional elec- 
tron gas (2DEG) (Ref. [J) marked evidence that incom- 
pressible states may form at unusual even-denominator 
filling fractions. After years of subsequent theoretical 
and experimental work^i^i^i^ it is well established that 
one of the most plausible theoretical candidates for a 
QH state at = | is an exotic state of matter, a paired 
quantum Hall state. Since electron-electron (e-e) inter- 
actions are repulsive this pair formation is a collective 
phenomenon involving residual interactions of composite 
particles that, in this state, are composites of an electron 
and two vortices. The electron pairing would be analo- 
gous to the formation of Cooper pairs in superconduc- 
tors, although it would be purely a result of e-e interac- 
tions without contribution from phonons or other fields. 
In some theoretical models, the excitations of the paired 
electron state are predicted to have non-abelian statistics 
that could be employed in the field of topological quan- 
tum computingi^ Currently, the most pressing challenge 
is to experimentally find evidence of the paired electron 
state and the particle statistics of its excitationsi^i^i^ The 
proposed testa^ for the non-abelian properties of these 
excitations make use of confined geometries and multi- 
ple constrictions in the 2DEG to generate interference 
among tunneling paths. This leads to a natural ques- 
tion whether the paired electron state is stable when the 
2DEG is confined into narrow trappings. 



This work addresses the structure of the v 



state 



when electrons in the 2DEG have been confined by exter- 
nal potentials into small QH droplets. They can be ex- 
perimentally realized by placing semiconductor quantum- 
dot (QD) devices into strong magnetic fields. "'^^ We show 
here theoretical evidence that in QH droplets the Pfaf- 
fian wave function,'^ which is commonly used to describe 



electron pairing, may have high overlaps with the exact 
many-body states at 1/ = |. In these calculations, we 
assume that the half-filled Landau level is spin-polarized 
and use realistic e-e potentials that include screening ef- 
fects from the background charge of electrons in the the 
lowest Landau level and a softening due to the finite 
thickness of the sample. However, the half-filled second 
Landau level of the Pfaffian state has a relatively high 
angular momentum, which may lead to its instability in 
the QD confinement. We show that in harmonic confin- 
ing potentials a compact filling of the half-filled Landau 
level is favored leading to the lowering of its angular mo- 
mentum. The paired electron state would then break 
down via fragmentation of spin and charge densities into 
two incompressible domains, spin-compensated i/ = 2 at 
the edges and spin-polarized = 3 at the center (see Refs. 
[TTI and . This phenomenon is analogous to the pro- 
posed formation of similar structures in the 2DEG where 
translational symmetry has been broken by long-range 
disorder. We present the fragmented states in QDs as 
alternatives to the Pfaffian state and show signatures of 
them in electron transport experiments. Based on these 
results, we conjecture that the stability of the paired elec- 
tron state depends crucially on the shape of the potential 
landscape where the electrons move in the 2DEG. This 
may explain, e.g., the observed fragility of the v — ^ QH 
state in narrow quantum point contactSf^^ 

The paper is organized as follows. We introduce our 
theoretical model of QDs in SecHlland the computational 
methods used to solve the many-body problem in Sec. lIIII 
The exact diagonalization method is used in Sec. IIVI to 
calculate the overlaps of the Pfaffian wave function with 
the exact many-body state. In Sec. |Vl we analyze the 
electronic structure of fragmented QH states and show 
that the second-lowest Landau level is spin-polarized due 
to the lifting of degeneracy of single-particle states near 
the Fermi level. In Sec. IVIl we present experimental ev- 
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idence for fragmentation of QH states in the 2 < < | 
filhng-factor regime. Section lVIII concludes our work with 
discussion of the relevance of our findings with the ob- 
served fragility of the v = ^ QH state in disordered or 
confined 2DEG. 



that both the DFT and QMC methods provide accurate 
results in this regime (see Appendix). 



A. Exact diagonalization 



II. MODEL 

QDs formed in the GaAs/Ala;Gai_a:As heterostructure 
are modeled for both lateral and vertical QD devices as 
droplets of electrons in a strictly two-dimensional (2D) 
plane confined by a parabolic external potential.— We 
use an effective-mass Hamiltonian 
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where N is the number of electrons, Vc{r) = rn'uj^r^ is 
the external parabolic confinement, m* — 0.067me is the 
effective mass and e — Yl.ltQ is the dielectric constant 
of the GaAs semiconductor medium, and A is the vector 
potential of the homogeneous magnetic field B perpen- 
dicular to the QD plane. The confinement strength ojq 
in the calculations is 2 meV, unless otherwise stated. 

If the e-e interactions are excluded, the single- 
particle solutions of the Hamiltonian ^ are Fock-Darwin 
stateSfiS, In the limit of a very high magnetic field, the 
Landau level structure approaches that of the 2DEG. 
However, in finite magnetic fields the external poten- 
tial alters the electronic structure and different Lan- 
dau levels overlap. Therefore, the concept of Landau 
level filling needs to be generalized to finite-size sys- 
tems. Kinaret ei al. defined the average filling factor 
as I'avg = N'^{N + L)/2, where L is the total angular 
momentum. Another possibility is to focus on the low- 
est Landau level (LLL) and define filling factor of a state 
as i^LLL — SA^/A'lll- These definitions differ in the high 
filling-factor regime, but this is not critical to the in- 
terpretation of results that are based on the structural 
properties of the many-body states. 



III. COMPUTATIONAL MANY-BODY 
METHODS 

The ground state corresponding to interacting elec- 
trons in QH droplets is solved numerically using the exact 
diagonalization (ED), density- functional theory (DFT), 
and the variational quantum Monte Carlo method 
(QMC). Since the paired electron state in the 2DEG is a 
strongly correlated many-body state, the ED method is 
used to analyze its stability in the QD confining poten- 
tial. The DFT and QMC methods are used to analyze 
the fragmented QH states. The regime where this frag- 
mentation gives characteristic signals in the experiments 
is beyond the reach of the ED method. However, we find 



In the ED method, we assume that the electrons oc- 
cupy states on one Landau level only. If we now take a 
fixed number of states from this Landau level, our com- 
putational task is first to construct the many-body basis. 
Then the Hamiltonian matrix corresponding to Hamilto- 
nian of Eq. IIJ is constructed in this basis. Finally, the 
lowest eigenstate and eigenvalue are found by matrix di- 
agonalization. More details can be found, e.g., in Ref . [TtI. 
In addition to the standard Coulomb interaction, we use 
in the ED two modifications of it. To model the finite 
thickness of the sample, we use a softened potentiali^ 
defined as 
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where dx is the sample thickness. Electrons in second or 
higher Landau levels move on top of background charge 
of lower Landau levels, which effectively screens the 
Coulomb interaction. This is modeled with a screened 
potential that is of the Gaussian form 
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where ds is the screening length. The unit of length 
in our ED results is given by I = ^ h/m*u)^ where 
Lo = ojq ~\- [ujc/i)"^ and = eB/m* is the cyclotron 
frequency of electron in magnetic field B. 



B. Density- functional theory 

Our DFT approach is based on spin-DFT, a variant of 
the conventional DFT generalized to deal with non-zero 
spin polarization. On top of standard spin-DFT, we in- 
clude the bare external vector potential A [see Eq. ([!])] 
in the Kohn-Sham equation. In contrast with current- 
spin-DFT, however, we neglect the exchange-correlation 
vector potential A^c- In the magnetic-field range con- 
sidered here, this has been shown to be a very reason- 
able approximationii^ As another valid approximation, 
we neglect the dependence of the exchange and corre- 
lation on the vorticityi^ The exchange and correlation 
energies and potentials are calculated using the 2D local 
spin-density approximation, for which we use the QMC 
parametrization of the correlation energy by Attaccalite 
et al?'^ 

The DFT approach is implemented on a 2D real- 
space grid and employs a multigrid method for solving of 
the Kohn-Sham equations. Our symmetry-unrestricted 
DFT approach has been shown to lead to solutions with 
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broken rotational symmetry that has been linked to mix- 
ing of the different eigenstates of angular momentumi^!.^'* 
In a fixed symmetric external potential, this type of spon- 
taneous symmetry breaking is expected to be unphysical. 
In Sec. IVIl we compare the validity of this assumption 
directly to experimental data. 

C. Quantum Monte Carlo 

Since the fragmentation of many-body state in the 
vicinity of = | is a delicate many-body problem, we 
employ the QMC method to analyze the reliability of our 
DFT approach. The wave function in the QMC is chosen 
to be 

N 

^ = D^D^l[j{n,) , (4) 

i<j 

where the two first factors are Slater determinants for the 
two spin types, and J is a Jastrow two-body correlation 
factor. The Slater determinants are constructed from the 
Fock-Darwin states. For the two-body Jastrow factor, we 
use a form 

^M=exp(^), (5) 

where a is fixed by the cusp condition to be three for 
a pair of equal spins and one for opposite ones, and 
b is an additional parameter different for both spin- 
pair possibilities. The ground state of the QD in the 
spin-droplet regime is calculated assuming that the LLL 
and the second-lowest Landau level (SLL) are com- 
pact. This means that the Slater determinants are 
built from single-particle states having angular momenta 
I = 0, . . . , A^LLL,s — 1 for spin s —1, | in the LLL and 
I = —1, . . . , iVsLL.s — 2 for the spin s in the SLL. En- 
ergy for each combination of non-negative total spin S 
and total angular momentum L is then calculated. The 
QMC method deals with the correlation effects in the 
many-particle system more accurately than the DFT ap- 
proach. However, the computational cost of the QMC 
is significantly larger than that of the DFT. A detailed 
description of the QMC method is given in Ref. [2^. 

IV. PFAFFIAN STATE IN QUANTUM DOTS 

The structure of the QH states in the 2DEG at half- 
integer filling factor has been a topic of intense research 
effortsi^ Currently, it is regarded plausible that the ex- 
perimentally observed = f state consists essentially 
of a full spin-compensated LLL and a half-filled spin- 
polarized SLL)^ in which weak p-wave electron pairing 
takes place. Formally, the SLL is described by a Moore- 
Read, or Pfafhan, wave function lifted to the SLLi^i^ 
There exists some theoretical evidence that the excita- 
tions of this QH state obey non-abelian statisticsi^i^ii^ 



ED calculations have become standard tests of trial wave 
functions of QH states, and they have shown high over- 
laps with the Pfafhan wave function in the 2DEGi^ 
However, there are other candidates for the v = ^ 
state, some of which possess only abelian quasiparticle 
excitations.— 

The structure of the v = \ state in QDs was analyzed 
with the ED method in Ref. illj. Here we provide results 
for half-filled higher Landau levels with more realistic 
inter-electron potentials defined in Sec. lIIII Following the 
theory of the v = ^ QH state in the 2DEG, we assume 
that the half-filled Landau level is spin-polarized. The 
Pfaffian wave function^S, which describes paired fermion 
states of the half-filled Landau level, is defined for LLL 
as 

In higher Landau levels the Pfaffian state is obtained 
by applying the Landau level raising operator to each 
electron. The angular momentum of the Pfaffian state is 
V = N'{N' - 1) - (nLL + ^)N', where N' is the number 
of electrons in the half-filled Landau level and tt-ll = 
{0, 1, . . .} is the Landau level index. 

We present the overlaps of the Pfafhan wave function 
with the ED eigenstate for electrons frozen to lowest 
(LLL), second (SLL), or third (TLL) Landau level, which 
correspond to filling fractions of = ^, v = |, and = |, 
respectively. Electrons in the half-filled second and third 
Landau levels move on top of the uniform background 
electron density of the spin-compensated lower Landau 
levels. This background charge effectively screens the 
Coulomb interaction. In QDs, the e-e interactions are 
further screened due to metallic leads. 

Figure [T] shows the overlaps of the Pfaffian wave func- 
tion and the ED eigenstate of Coulomb interaction for 
particle numbers 4 < TV' < 12. For large particle num- 
bers, the overlaps in the second Landau level are highest. 
This shows that — ^ has the highest probability to be 
described by the Pfafhan. 

Next, we study how the screening of the e-e interaction 
and finite thickness of the sample change the overlaps of 
the ED eigenstate with the Pfafhan. For six electrons on 
LLL, the overlaps are slightly improved when the screen- 
ing and finite sample thickness are taken into account in 
the interaction (see Fig. [21). On SLL, screening slightly 
improves the overlap, but a finite thickness lowers it. The 
same trends can be seen in Fig. [5] for eight electrons, but 
now the effects are clearly stronger, and there is a large 
increase in the overlaps. On the LLL, a finite sample 
thickness is needed to achieve the best overlap. On the 
SLL, the screening increases the overlap, which can be 
contrasted with the spherical geometry where the SLL 
overlap is maximized at a finite thickness of the samplei^S 

The highest overlaps are of the order of 0.8-0.9 at 
!/ = |, which means that the structure of the many- 
body state is close to the Pfaffian. The exact state at the 
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FIG. 1: (Color online) The overlaps of the Pfaffian wave func- 
tion with the corresponding exact state at the lowest (LLL), 
second (SLL), and third (TLL) Landau level in the case of 
a Coulombic electron-electron interaction. A^' denotes the 
number of electrons in the half-filled Landau level. 

given angular momentum would therefore show electron 
pairing to a significant degree. We note that the Pfaffian 
wave function in Eq. ^ has no adjustable parameters. It 
is possible to modify the Pfaffian wave function by intro- 
ducing a pairing function that differs from g = — Zj ) 
of the Moore- Read formf22, This won't change the angular 
momentum of the state but has been found to increase 
overlaps significantly in the 2DEG. 

In addition to the overlaps, it is crucial to verify that 
the state at the angular momentum of the Pfaffian state 
is energetically favorable. In fact, the LLL = | state 
with A^' = 8 corresponding to Fig. ^d) is a possible 
ground state at small and large values of the thickness 
dx, but not at values of dx where the overlap is peaked j2i 
A further obstacle for the Pfaffian state in finite-size QH 
droplets is that the SLL may not attain the high angular 
momentum and complete spin polarization of the Pfaf- 
fian. In QH droplets, the degeneracy of Landau levels is 
lifted when electrons move in external confining poten- 
tials [Fig. [SJb)], and a compact distribution of electrons 
on the Landau levels could be energetically more favor- 
able. In the next section, we show that this would lead 
to nonexistence of the paired electron state and intro- 
duce fragmented QH states in quantum Hall droplets as 
alternatives. 



V. FRAGMENTED QUANTUM HALL STATES 

In quantum Hall droplets, single-particle states within 
each Landau level are not degenerate due to the confin- 
ing potential. The average distance of an electron from 
the center of the droplet, and therefore also the potential 
energy, increase with angular momentum. This suggests 



that a compact occupation structure may be energeti- 
cally favorable. The compact occupation of Landau lev- 
els leads to fragmentation of charge and spin densities 
into incompressible integer filling factor domains. We 
call these states fragmented quantum Hall states that are 
alternatives to the paired electron state at half-integer 
Landau level fillings. 

We analyze the structure of fragmented QH states near 
J/ = I in a harmonic confining potential of a semicon- 
ductor quantum dot with the QMC and the DFT meth- 
ods. The Kohn-Sham single-particle energy spectrum of 
the Landau levels calculated with the spin-compensated 
DFT and the spin-DFT are shown in Fig. [3)^a) and (b), 
respectively. The spin-DFT and the QMC show that the 
degeneracy of the single-particle states close the Fermi 
energy is lifted via a complete polarization of the second- 
lowest Landau level. Therefore, a compact occupation 
of the single-particle states of the spin-compensated LLL 
and spin-polarized SLL leads to a fragmented state with 
a. v — 2 region (double-occupied LLL) at the edges of 
the droplet and v = 3 (spin-polarized SLL) at the center 

[Fig.m. 

The spin-splitting of the SLL in the spin-DFT calcu- 
lations is analogous to the Stoner criterion, which states 
that in the presence of correlations between electrons of 
the same spin and high density of states near the Fermi 
level, the system prefers ferromagnetic alignment that 
reduces the degeneracyi^S We call the incompressible, 
spin-polarized droplet of SLL electrons a spin droplet— 
The size of the spin-droplet becomes significant when the 
number of electrons in the dot N > 35. The non- uniform 
filling-factor structure of the spin-droplet states is rem- 
iniscent of the incompressible QH domains that form in 
the 2DEG with long-range disorderJ^ The compact oc- 
cupation of the SLL leads to a lower angular momentum 
than what is needed for a paired electron state as de- 
scribed by the Pfaffian wave function For example, 
the size of the spin-droplet in the QMC method is TV' = 8 
electrons at A'' = 48, and the angular momentum of the 
SLL is L' = 20, which can be contrasted to L' — 44 
for the Pfaffian wave function with the same number of 
electrons. 

The SLL remains polarized and compact between | > 
1/ > 2. Hence, we call this filling-factor range the spin- 
droplet regime. The size of the spin-droplet gradually 
shrinks with the increasing magnetic field as the electrons 
are passed from the SLL to the LLL. The contributions 
of the LLL and SLL occupancies to the electron and spin 
densities are shown in Fig. [4] for the case of 60-electron 
QD. Qualitatively similar results were obtained for con- 
finement strengths 1 to 4 meV and electron numbers N 
between 35 and 120, which conffims the generality of the 
results. The calculations show that the energy benefit 
from the polarization of the SLL is large (see Fig. [3] and 
Appendix), which would make spin-droplets robust in the 
presence of impurities in samples. 

We note that the stability of the fragmented QH states 
in large QDs {N > 30) can be contrasted to the insta- 
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FIG. 2: Overlap of the Pfaffian wave function with the corresponding exact state for TV' = 6 and N' = 8 electrons in the half- 
filled Landau level, respectively, using (a-b) screened electron-electron potential with screening length ds and (c-d) softened 
potential due to finite sample thickness dr for electrons at the lowest Landau level (LLL) corresponding to ly — i, and the 
second lowest Landau level (SLL) corresponding to = | . 



bility of the maximum-density-droplet (MDD) state in 
the same regime. The MDD state is the totally polar- 
ized state corresponding to the 1^=1 QH state in 2DEG, 
and it has been found to be instable in large QDs with 
iV > 30 (Refs. [Hand [11). 



VI. SIGNATURES OF FRAGMENTATION IN 
ELECTRON AND SPIN TRANSPORT 

The emergence of finite size counterparts of integer and 
fractional QH states in QDs gives characteristic signa- 
tures in the chemical potentials. Several experimental 
methods have been developed to measure the chemical 
potential in a QD via addition of electrons one-by-one 
into the system. These experimental methods include 
Coulomb blockade;^ capacitance?^ and charge detec- 
tion techniquesi^ In this work, we use data from elec- 
tron transport measurements of QDs in the Coulomb and 
spin blockade regimei^i^^ The spacings of the spin and 
Coulomb blockade peaks correspond to the energy needed 
to add the A^th electron in the system of iV — 1 elec- 
trons, i.e., the chemical potential defined as fi{N,B) — 
EiotiN,B)-Etot{N-l,B). 

We calculate the signatures in the chemical potentials 
associated with the formation of fragmented QH states 



and compare these to those obtained from the electron 
transport data in three different QD devices. Two of 
the experimental samples (sample A and B) are lateral 
quantum dots on a high-mobility 2DEG^ while the third 
one (sample C) is a vertical QD.'^^ The samples A and 
B were manufactured on a high mobility 2DEG samples 
with spin-polarized leads for electron transport measure- 
ments in the spin blockade regime. The data of the sam- 
ple C was obtained in the Coulomb blockade. The high 
mobility of samples chosen for comparison is essential to 
reduce unpredictable effects of impurities and disorder 
that make identification of signals of physical phenom- 
ena difficult. 

We first address the problem of whether the electronic 
states in the QD samples show any signs of broken rota- 
tional symmetry. Inhomogeneities and impurities in QD 
devices may break the rotation symmetry, and a Jahn- 
Teller type of mechanism could be active if disorder alters 
significantly the shape of the confining potential. As a 
result, the ground state transitions with increasing mag- 
netic field become continuous rather than discrete. A 
signature of this type of symmetry breaking would be a 
smoothing of the chemical potential. Experimental data 
from a high-mobility lateral QD device is of sufficiently 
good quality to test for the presence of symmetry break- 
ing mechanisms. Figure [5] shows a comparison of the elec- 
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FIG. 3: (Color online) (a) Kohn-Sham energy spectrum of a 
60-electron quantum dot as a function of single-particle angu- 
lar momentum L calculated from the density functional the- 
ory with spin-compensated orbitals. The density of states of 
the second-lowest Landau level (SLL) is high near the Fermi 
energy _Ef. The magnetic field is 2.125 T which corresponds 
to a filling factor oi v = |. (b) The corresponding energy 
spectrum from spin-density functional theory shows lifting of 
the degeneracy near the Fermi level via complete spin polar- 
ization of the SLL. The lowest Landau level (LLL) remains 
spin compensated. Spin | (J.) corresponds to spin orientation 
parallel (antiparallel) to the magnetic field. The spin-splitting 
due to many-body effects is about 1.5 meV at L = 0. In com- 
parison, the Zeeman splitting is about 0.05 meV. 
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FIG. 5: (Color online) Ground state energy of the N = Al 
quantum dot calculated with the density-functional theory for 
the rotationally symmetric eigenstates of the angular momen- 
tum L (blue curve) and for the ground states in the symmetry- 
unrestricted approach (magenta curve). The corresponding 
chemical potentials fj,{N = 47 — > 48) calculated from the the- 
ory are shown in the lower panel together with experimental 
data from a lateral quantum-dot device (sample B). Dashed 
lines correspond to the boundaries of the spin-droplet regime. 
The insets show the fragmented spin and charge densities of 
three of the corresponding states (cf. Fig. The strength 
of the parabolic confining potential of the quantum dot is 
ijjo = "2 meV in the calculations. 



FIG. 4: (Color online) Total electron density p-\ -j- pi (full 
region) and the net spin density p| — pi (transparent blue 
region) of quantum Hall states in a quantum dot at (a) /^ = 2, 
at (b) an intermediate state between u = 2 and = |, and at 
(c) V = |. The latter two show fragmented charge and spin 
densities with spin-compensated v = 2 region at the edges 
and spin-polarized = 3 at the center. The densities were 
calculated with the spin-density-functional theory for a 60- 
electron quantum dot. The net spin-up density is due to spin 
polarization of the second-lowest Landau level. 



tron transport data to DFT calculations with and with- 
out symmetry breaking. The data shows sharp increases 
in the chemical potentials, which is consistent with dis- 
crete transitions in the ground state. Therefore, to a good 
approximation, the rotational symmetry is preserved in 
high mobility samples, and the angular momentum L is 
a good quantum number. 

The complete polarization of the SLL at = | is re- 
flected in the energetics of the system. The DFT cal- 
culations show that at = I there is a step feature 
followed by a plateau region in the chemical potential. 



Figure [5] shows the DFT results for chemical potentials 
of = 24 . . . 48 in comparison with the experiments. 
The step in the chemical potential is associated with the 
total polarization of the SLL in the DFT calculations. 
This feature can be found in all three samples we studied 
above A'" > 30, as predicted by the theory. Some models 
of QDs assume that the SLL is spin-polarized due to the 
Zeeman effect. This model does not, however, apply for 
the lateral and vertical QD devices examined in this work 
where the effect of the Zeeman splitting is estimated to 
be only a few percent of the spin splitting caused by the 
many-body interactions (see Fig. 

In the ^ > u > 2 regime, the ground state energy is 
approximately constant (see Appendix), and the calcu- 
lations show a phase transition in the system where two 
phases (i/ = 2 and v = ?>) coexist, and the size of the 
V = 2 domain increases with the magnetic field. The 
chemical potential does not continue to rise, but instead, 
it is oscillating around a constant value until v = 2. 
This signature can be found in all the experimental sam- 
ples (see Fig. [5]). All electron transport data presented 
are thus consistent with the theoretical picture that the 
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FIG. 6: (Color online) Chemical potentials calculated with the density-functional theory (a) and measured from a vertical 
(b) and lateral (c) quantum-dot devices for various electron numbers. Both experiments show the signal associated with the 
polarization of the second lowest Landau level at i/ = 5/2 in the peak position data when A'^ > 30 in agreement with the 
theoretical result. The confinement strength ranged from 2 to 4 meV depending on the electron number. The data for the 
vertical device in (b) is courtesy of L. Kouwenhoven,— and the data for lateral device in (c) is courtesy of A. S. Sachrajda.— 



ground states in the vicinity of v — ^ involve fragmented 
QH states. We point out, however, that the results are 
sensitive to the shape of the external potential, and the 
pairing of the electrons may still occur if the potential is 
sufficiently homogeneous, e.g., in large QDs, where the 
second Landau level would acquire higher angular mo- 
mentum. 

Spin polarization of the leads is commonly used to cre- 
ate a current that depends on the orientation of the elec- 
tron spin, which passes through the device. In the case 
of the two lateral QDs in our analysis, the electrons en- 
ter the QD from spin-polarized magnetic edge states of 
the 2DEG through tunneling barriers. Coulomb block- 
ade lifts when the energies of the many-body states cor- 
responding to N and -I- 1 electrons are equal. The 
tunneling current depends then on the coupling between 
the wave function in the QD and the electronic states 
in the external leads. The lowest Landau level orbitals 
are at the edges of the QD, and the coupling is stronger 
to the leads compared to the second lowest Landau level 
orbitals that are close to the center of the QD. Due to 
polarization of the leads, their coupling to electron states 
with spin parallel to the external polarization is higher 
than the coupling of spins antiparallel to the external 
polarization. This spin dependence in the transport has 
been shown to lead to a characteristic checkerboard pat- 
tern of current densities through QDsi^'^'^'^i^ Our 
DFT results are consistent with such transport currents 
in the spin blockade regime (Fig. [7]). The polarization of 
the SLL in the | > > 2 regime would be in contrast to 
the model presented in Ref. Sol. A consequence of this is 



that the transport current via SLL orbitals should show 
no checkerboard pattern in this regime since the spins 
are always parallel to the polarization of the leads. This 
could be tested with high-accuracy spin blockade spec- 
troscopy which would be able to detect small changes in 
the weak tunneling currents through the SLL orbitals. 



VII. FRAGILITY OF THE = | QUANTUM 
HALL STATE 

The V — ^ state is one of the most fragile QH states. 
It is observed only in high-mobility 2DEG samples as the 
paired electron state may break down in the presence of 
impurities. These induce a non-uniform potential that, in 
the light of results in this work, may lead to its instabil- 
ity. Our findings are thus in line with those obtained by 
Chklovskii and Lee who predicted that in the presence of 
long-range disorder in the 2DEG, incompressible integer 
filling factor regions form that are separated by domain 
walls These structures are analogous to the fragmented 
QH states that we find in QDs. Structures reminiscent 
of domain walls have been observed with scanning-probe 
imaging techniques in a perturbed QH liquid4^ 

Analogous instability of QH states may also occur in 
other geometries where the electrons are not strictly con- 
fined in all directions, such as in high-mobility 2DEG 
samples in the vicinity of constrictions. One indication 
of this may be the observed fragility of the v = ^ state 
in narrow quantum point contacts.^* Proposed tests^ for 
the non-abelian properties of quasi-particle excitations of 
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DFT calculation 



N=48 



o-a N=47 / 



Experiments (Sample A) 
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N=46 
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FIG. 7: (Color online) Checkerboard pattern of transport current in density-functional theory (left panel) and spin-blockade 
experiments (right panel) . The lowest current densities correspond to electron transport via states in the second-lowest Landau 
level, near the core of the quantum dot. The current density in experiments has been amplified in high magnetic fields with a 
linear function to compensate for the general attenuation of the signal. 



— J QH state make use of finite geometries and mul- 
tiple constrictions to generate interference among tun- 
neling paths. A possible fragmentation of the v = ^ 
QH state close to the boundaries, which would lead to 
the instability in such geometries, is still an open ques- 
tion that requires further analysis of the effects of the 
confinement. While recent experiments on the quasi- 
particle tunneling,^ shot noise generated by partitioning 
edge currents,— and interferometric measurements of QH 
edge excitations^ of the v = ^ QH state show results, 
which are all consistent with the unusual quasi-particle 
charge e* = | of the paired electron state, the parti- 
cle statistics of the excitations remains to be confirmed. 
Possible fragmentation of QH states in narrow constric- 
tions needed for quasi-particle interferometry adds an- 
other challenge in this long quest to confirm the possible 
non-abelian characteristics of the = | state. 

To conclude, we have shown theoretical evidence that 
electron pairing is possible in small QH droplets in quan- 
tum dots at = |, provided that the half- filled Lan- 
dau level can acquire sufficiently high angular momen- 
tum. However, our calculations indicate that in parabolic 
external confining potentials the paired electron state 
breaks down leading to fragmentation of charge and spin 
densities. We find indirect evidence of such fragmenta- 
tion in several experiments but point out that our results 
can be tested by direct measurements of the spatial de- 
pendence of spin and charge densities in different geome- 
tries and experimental setups. 
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VIII. APPENDIX: ACCURACY OF 
NUMERICAL METHODS 

The electron correlations play an important role in the 
structure of fractional QH states. To test for the accu- 
racy of the DFT method in the spin-droplet regime, we 
compare the energies of different spin polarization states 
between the DFT and the QMC in the | > > 2 regime. 
The results for a 48-electron QD are shown in Fig. [S] 

Both methods show the spin-droplet structure with a 
comparable energy benefit in the polarization (5 « 0.5 
meV for S'max = 4. The QMC method estimates that the 
maximum size of the spin droplet is iVsD — 7 compared to 
8 in the DFT. Given the typical statistical error of ±0.05 
meV in the QMC results, the overall agreement between 
the methods is excellent. This test indicates that the 
DFT method captures the essential many-body physics of 
the spin-droplet formation and gives accurate results for 
the ground states. The DFT method was subsequently 
used in the calculation of the chemical potentials of large 
QDs, which can be compared to the transport experi- 
ments in the spin blockade regime. 

The DFT method predicts some non-compact states 
outside the spin-droplet regime, e.g., L — 375, S — 7/2 
state as shown in Fig. [5] This state has one spin-down 
electron in the SLL with Z = 0. Emergence of non- 
compact states is a manifestation of the degeneracy of 
the single-particle states near Fermi-level. However, they 
are rare in the DFT and occur only at magnetic fields 
below the polarization of the SLL. Detailed analysis of 
these states with the QMC goes beyond the scope of the 
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present work and is left for future research. 



1 1.2 1.4 1.6 

B(T) 

FIG. 8: (Color online) Comparison of the ground state en- 
ergies for given total spin S in the density-functional theory 
(DFT) and the quantum Monte Carlo (QMC) method. The 
number of electrons A'^ = 48. The line widths in the QMC 
denote the statistical error in the results. Only the ground 
state and the S = state are shown in the DFT result. The 
strength of the parabolic confining potential of the quantum 
dot is Wo = 2 meV in the calculations. 
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